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Abstract. This paper reports on the methods and results of a theoretical 
analysis to design an insulator which must provide a thermally quiet environment 
to test on ground delicate temperature sensors and associated electronics. These 
will fly on board ESA's LISA Pathfinder (LPF) mission as part of the thermal 
diagnostics subsystem of the LISA Test-flight Package (LTP). We evaluate the 
heat transfer function (in frequency domain) of a central body of good thermal 
conductivity surrounded by a layer of a very poorly conducting substrate. This 
is applied to assess the materials and dimensions necessary to meet temperature 
stability requirements in the metal core, where sensors will be implanted for test. 
The analysis is extended to evaluate the losses caused by heat leakage through 
connecting wires, linking the sensors with the electronics in a box outside the 
insulator. The results indicate that, in spite of the very demanding stability 
conditions, a sphere of outer diameter of the order one metre is sufficient. 
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1. Introduction 

LISA Pathfinder (LPF) is an ESA mission, whose main objective is to put to test 
critical parts of LISA (Laser Interferometer Space Antenna), the first space borne 
gravitational wave (GW) observatory pp. The science module on board LPF is the 
LISA Test- flight Package (LTP) which basically consists in two test masses in 
nominally perfect geodesic motion (free fall), and a laser metrology system, which 
reports on residual deviations of the test masses' actual motion from the ideal free 
fall, to a given level of accuracy [jj]. 

In order to ensure that the test masses are not deviated from their geodesic 
trajectories by external (non-gravitational) agents, a so called Gravitational Reference 
System (GRS) is used This consists in position sensors for the masses which send 
signals to a set of micro-thrusters; the latter take care of correcting as necessary the 
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spacecraft trajectory, so that at least one of the test masses remains centred relative to 
the spacecraft at all times. The combination of the GRS plus the actuators is known 
as drag- free subsystem||. 

The drag-free is of course a central component of LISA, and needs to be operated 
at extremely demanding levels of accuracy. The laser metrology system should then 
be sufficiently precise to measure relative test mass deviations. The overall level of 
noise acceptable for LISA is defined in terms of rms acceleration spectral density, and 
has been set to 
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in the frequency range fO 4 Hz < uj/2it < 10 3 Hz. This is equivalent to 
Sj/ 2 ~ 4xI0~ 21 Hz -1 / 2 , with the same frequency dependence. 

Because LPF is a technological mission, aimed to assess the feasibility of LISA, 
its ultimate goal has been relaxed to [Ej 
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in the frequency range 1 mHz < lu/2tt < 30 mHz, i.e., one order of magnitude less 
demanding, both in noise amplitude and in frequency band. 

Equation J5J gives the global noise budget. This is naturally made up of 
contributions from different perturbative agents, such as temperature and magnetic 
field fluctuations, GRS and interferometer noise, etc. As a general rule, a requirement 
on the magnitude of each of the various perturbing factors is set at a I % fraction of 
the total. In the case of temperature fluctuations, this is equivalent to 

4 /2 H< ir4K /^. lmHz < w/2tt < 30mHz (3) 

Because temperature stability is important, a decision has been taken to place 
high precision thermometers in several strategic spots across the LTP — as part of what 
is called Diagnostics Subsystem [S| %■ Such high precision temperature measurements 
will be useful to identify the fraction of the total system noise which is due to thermal 
fluctuations only, and this will in turn provide important debugging information to 
assess the performance of the LTP. 



1.1. Temperature measurements 

If the temperature gauges are to be sensitive to fluctuations at the level given by 
P| then clearly the entire measuring device should be less noisy, typically by a factor 
of 10. This means that such device, which includes both the sensors and the associated 
electronics, can generate a maximum level of noise of 

5 7> e „sor M < lcr5 K/VHz , I mHz < uj/2it < 30 mHz (4) 

Research work is currently being conducted at IEEC (Barcelona, Spain) to 
identify the appropriate sensors and design the better suited front end electronics. But 

| The term drag-free dates back to the early days of space navigation, when it was used to name a 
trajectory correction system designed to compensate for the effect of atmospheric drag on satellites 
in low altitude orbits. 

% The Diagnostics Subsystem of the LTP also includes magnetometric measurements and a charged 
particle flux detector. 
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the prototype system needs of course to be tested for compliance with equation ijljl. 
Thus, in order to do a meaningful test, the system must be sufficiently thermally 
isolated that the observed fluctuations in the readout data can be attributed solely 
to sensor noise, rather than to a combination of it with real ambient temperature 
fluctuations. This means temperature fluctiations in the thermomters' placements 
should again be at least one order of magnitude below the target sensitivity, 
equation Q, or 

S^ esthed (u) < I CT 6 K/Vlz , f mHz <lo/2tt< 30 mHz (5) 

It turns out that 10~ 6 K/v / Hz is a truly demanding temperature stability, orders 
of magnitude beyond the capabilities of normal thermally regulated rooms. We 
thus need to design a specific thermal insulator to shield the sensors from ambient 
temperature fluctuations during the test process. 

In the ensuing pages we describe in detail the insulator design. It is extremely 
important to stress at this point that the performance of the insulator, i.e., its ability 
to screen out ambient temperature fluctuations, cannot be checked experimentally, 
at least under working thermal conditions in the laboratory. This is because, by 
definition, the insulator is the tool to check the sensing instruments, not viceversa: we 
need to rely on the results of theoretical argumentation to make a decission on which 
is the appropriate thermal insulator for our purposes. An experimental verification 
of the model is only thinkable under much more extreme conditions, where external 
temperature fluctuations are orders of magnitude higher than the ones which will be 
met during the test. 

2. Thermal insulator design concept 

The idea of the insulator design is displayed in figure^ an interior metal core of good 
thermal conductivity is surrounded by a thick layer of a poorly conductive material. 
The inner block ensures thermal stability of the sensors attached to it, while the 
surrounding substrate efficiently shields it from external temperature fluctuations in 
the laboratory ambient. We propose a spherical shape for the sake of simplicity of the 
mathematical analysis, even though this will be eventually changed to cubic in the 
actual experimental device due to practical feasibility issues. 

2.1. Mathematical model 

The basic assumption of the mathematical analysis we shall present is that heat flows 
from the interior of the insulator to the air outside, and from the latter to the interior 
of the insulator, only by thermal conduction. This is a very realistic hypothesis in 
the context of the experiment, as radiation mechanisms are certainly negligible and 
convection should not play any significant role, either, since the entire body is solid, 
and temperature fluctuations will be small at all times anyway. 

Let then T(x, t) be the temperature at time t of a point positioned at vector x 
relative to the centre of the sphere. T(x, t) thus satisfies Fourier's partial differential 
equation [7] 

pc p ^T(x,i)=V-[ K VT(x,i)] (6) 

where p, c p and k are the density, specific heat and thermal conductivity, respectively, 
of the substrate. We shall assume these are uniform values within each of the two 
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materials making up the insulating body, with abrupt changes in the interface. We 
can thus represent them as discontinuous functions of the radial coordinate, as follows: 

I P2,C P 2,K2 It 0-1 S T < a-2, 

with r = |x|. Initial and boundary conditions are the following: 

T(x,i = 0) = 0, T(r = a 2 ,t) = T (6,ip;t) (8) 

where 9 and tp are spherical angles which define positions on the sphere's surface. The 
boundary temperature can be expediently expressed as a multipole expansion: 

T {9, cp; t) = hm (t) Y lm (9, ip) (9) 

1 m 

where Yi m (0,ip) are spherical harmonics, and bi m (t) are boundary multipole 
temperature components. 

In practice, the boundary temperature will be randomly fluctuating, therefore 
bi m (t) will be considered stochastic functions of time. We shall also reasonably assume 
them to be stationary Gaussian noise processes with known spectral densities, Si m (uj). 

As shown in the appendix, the frequency analysis of this problem leads to a 
transfer function expression of the temperature inside the body: 

T(x, lj) = ^ Hlm(x, bi m (u) (10) 

I in 

where tildes (~) stand for Fourier transforms, e.g., 

T(x,uj)= T(x,t)e- iiAjt dt (11) 



etc. If we make the further assumption that different multipole temperature 
fluctuations at the boundary are uncorrelated, i.e., 

(b* lm ,{uj) b~im(u)) = Si m (aj) Sin S m > m (12) 




Figure 1. Thermal insulator design concept. Left: 3D diagramme, including 
sensor placement principle. Right: cut view, with notation convention dictionary. 
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then the spectral density of fluctuations at any given point inside the insulating body 
is given by 



St(x,w) = \Hi m (x,Lu)\ 2 Si m (u>) 



(13) 



It is ultimately the spectral density St(x-, uj) which has to comply with the 
requirement expressed by equation Q). Based on knowledge (by direct measurement) 
of ambient laboratory temperature fluctuations, equation l|I3|l will provide the 
guidelines, as regards materials and dimensions, for the actual design of a suitable 
insulator jig. 



3. Homogeneous boundary conditions 

Thermal conditions in the laboratory are rather homogeneous. This means that the 
boundary temperature fluctuations will be in practice essentially independent of the 
angles and ip, i.e., 

T (fi,wt)=B{t) (14) 

and consequently the generic expansion equation 10 includes only the monopole term, 
hence 



b 00 (t)=V4TrB(t) (15) 

The temperature T(x, uj) in this case will only depend on radial depth, r, 
therefore, 

f{r,u>) = H(r,u)B(u>) (16) 
with H (r, uj) = -\/47r i?oo( x , ^)- According to equation (|A. 2 1|> of the Appendix, this is 



H(r,u) = 



Vo(u) joi^r) + Co M ?/o (72?") j a i < r < a 2 



(17) 
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Figure 2. Frequency response at the centre (r = 0) of a spherical thermal 
insulator — as given in equation 1171 — , along with the frequency response of 
a low-pass filter of order one and same frequency cut-off, u; cut _ ft = 1/r. 
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This is a low-pass filter transfer function — even though the cumbersome 
frequency dependencies involved in the expressions above do not make it immediately 
obvious. A plot of the square modulus of H(r,u>) is shown in figure |21 for r = (red 
curve) . The figure also shows a low-pass filter of the first order with the same frequency 
cut-off, |-ffist ordcr(w)| 2 = (1 +tj 2 T 2 )~ 1 , for conceptual comparison (blue curve). 

The most salient feature emerging out of the plot is the stronger drop in H(r, uj) 
at the high frequency tails. The latter can be easily assessed in quantitative detail, 
and the result is 

\H(Q,u)\ ^wre"^ 7 (18) 

where r is the filter's time constant — a complicated function of the insulator's physical 
and geometric properties, to be discussed below. 

As already mentioned in the Introduction section, to test the temperature sensors 
and electronics we need a very strong noise suppression factor in the LTP frequency 
band. A look at figure |21 readily shows that high damping factors require such 
frequency band to lie in the filter's tails. The thermal insulator should therefore 
be designed in such a way that its time constant r be sufficiently large to ensure 
that the LTP frequencies are high enough compared to 1/r. The exponential drop in 
the transfer function shown by equation (|18|l makes the filter actually feasible with 
reasonable dimensions. 

4. Numerical analysis 

In this section we consider the application of the above formalism to obtain practically 
useful numbers for the actual implementation of a real insulator device which complies 
with the needs of our experiment. 

First of all, a selection of an aluminum core surrounded by a layer of polyurethane 
was made. Aluminum is a good heat conductor and is easy to work with in the 
laboratory; polyurethane is a good insulator and is also convenient to handle, as 
it can be foamed to any desired shape from canned liquid. Other alternatives are 
certainly possible, but this appears sufficiently good and we shall therefore only make 
reference to this specific one. 

The relevant physical properties of aluminium and polyurethane are specified in 
table □ 

Density Specific heat Thermal conductivity 
p(kgm- 3 ) c p (Jkg^K" 1 ) ^(Wm-iK- 1 ) 



Aluminum 2700 900 250 

Polyurethane 35 1000 0.04 

Table 1. Density, specific heat and thermal conductivity of aluminium and 
polyurethane. Units are given in the International System. 

Figure plots the amplitude damping coefficient of the insulator block, \H(r,u>)\, 
at the lower end of the LTP frequency band, i.e., 1 mHz, and at the interface position, 
r = a%. Each of the curves corresponds to a fixed value of the latter, and is represented 
as a function of the outer radius of the insulator. This choice is useful because the 
sensors are implanted for test on the surface of the aluminium core, and also because 
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at higher frequencies thermal damping is stronger. So in practice the actual damping 
power of the device will be the one plotted, and better at the higher frequencies in 
the measuring bandwidth. The figure clearly shows that the assymptotic regime of 
equation (|18|l is quite early established. 

The choice of dimensions for the insulating body must of course ensure that the 
minimum requirement, equation (|SJ) is met. For this, a primary consideration is the 
size of the ambient temperature fluctuations in the site where the experiment is made. 
Dedicated measurements in our laboratory showed that 

S^LbicntH - 10" 1 K/Vlz , 1 mHz < w/2tt < 30 mHz (19) 

We therefore need to implement a device such that \H{cl\,uj) \ < 10~ 5 throughout 
the measuring bandwidth (MBW). Suitable dimensions can then be readily read off 
figure |3 and various alternatives are possible, as seen. Before making a decission, 
however, we need to make an additional estimate of the heat leakage down the electric 
wires which connect the temperature sensors with the elctronics, which lies of course 
outside the insulator. We come to this next. 

4-1. Heat leakage through connecting wires 

We use a simple model, consisting in assuming the connecting wires behave as straight 
metallic rods which connect the central aluminum core with the electronics, placed 
in the external laboratory ambient. Because the polyurethane provides a very stable 
insulation, we can neglect the lateral flux, hence only a unidimensional heat flow 
needs to be considered. For this, the following equation relates the heat flux to the 
temperature difference between the two wires' edges: 

ttR 2 - 

Q(t) = Kwirc -j^ [T(a 2 ,t) - T(ai,t)] (20) 

^-wirc 

where K w ; re is the thermal conductivity of the wire, i? W iro its transverse radius, and 
£ w i rc its length inside the polyurethane layer. 
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On the other hand, the heat flux results in temperature variations in the metal 
core, given by 

dT 

Q(t) = Pic pl V 1 —(a 1 ,t) (21) 

where V\ = Aira\/'S is the volume of the metal core. Equating the above expressions 
we find 

K w i Ie ^^[T(a 2 ,t)-T(a 1 ,t)]=p 1 c pl V l ^ r (a u t) (22) 

For fluctuating temperatures, we can now obtain the relationship between the 
spectral density at the aluminium core and the ambient, due to heat conduction along 
the wire: 

■^wircK' W ) = l^wircHI ^ambient H ( 23 ) 

where 

\H wile (")\-- * wire T ff e (24) 

and where the approximation has been made that the temperature fluctuations at the 
inner end of the wire are much smaller than those at the outer end, due to the presence 
of the polyurethane layer. 

In practice, there will be several sensors for test inside the insulator. Under the 
hypothesis made that no lateral heat flux is relevant, the transfer function for a bundle 
of N of wires is, at most, N times that of a single wire. Thus, 

Iff I \\ I — ^ ^wirc -Rwire ( r\ 

-H JV wires M - — - — — (25) 

W/27T 87TjOiCpiaf «wiro 

Let us consider numerical values in this expression. We use thin copper wires 
(kcu = 401 Wm _1 K _1 ) of radius i? w i re = 0.1 mm, and assume some fiducial parameters 
for the size of the aluminium core, a\, the wire length, £ w i re , the number of connecting 
wires, N, and the frequency, lj/2tt. The following obtains: 



\H NwiI M\ = i.i x io- 5 f^V^rVi^Vf— V 1 <*» 

1 Nwncsy n V 3 °/^ 13cm ^ v 25cn v v lmHz / 

This result indicates that, for laboratory fluctuations in the level of equation 119|) . 
leakage through wiring causes fluctuations in the temperature sensors of about 
10 _6 K/v / Hz, equation (J2HJ), which is compliant with the requirement of stability of 
equation JSJ. The most sensitive parameter in the above expression is the size of the 
metal core, and this determines the need to make it somewhat large. The length of 
the wires has been taken to be 25 cm, but this does not necessarily mean we need 
<i2 = 38 cm (assuming the radius of the aluminum core is a\ = 13 cm), because the 
wires can be partly wound inside the polyurethane layer to further protect the system 
against leakage. In fact, this wire lengthening is an easy way to improve attenuation. 

As regards frequency dependence, compliance is guaranteed in the entire MBW if 
it is at its lower end: indeed, not only |iT W j re (a;)| decreases as also ambient noise 
fluctuations drop below 10 _1 K/VHz at higher frequencies. 
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5. Conclusions 

Temperature fluctuation measurement is very demanding in the LTP, and 
subsequently LISA, as reflected by equation J3J|. Accordingly, very delicate sensor 
and associated electronics must be designed, and of course tested in ground before 
boarding. 

However, even the best laboratory conditions are orders of magnitude worse 
than the above requirement, so meaningful tests of the temperature sensing system 
cannot be tested without suitably screening the sensors from ambient temperature 
fluctuations. We have addressed how this can be accomplished by means of an 
insulating system consisting of a central metallic core surrounded by a thick layer 
of a very poorly conducting material. The latter provides good thermal insulation, 
while the central core, having a large thermal inertia, ensures stability of the sensors' 
environemnt. The choice of materials is flexible, so aluminium and polyurethane, 
which are easily available in the market, has been adopted. Thereafter, the dimensions 
need to be fixed. 

The appropriate sensors for the needs are temperature sensitive resistors, more 
specifically thermistors — also known as NTCs. It appears that, because these sensors 
need to be wired to external electronics, heat leakage through such wires is an effect 
which needs to be quantitatively assessed to prevent losses. We have analysed this 
problem, and concluded that it strongly depends on the central metallic core size, and 
imposes that it be somewhat large. 

Laboratory ambient temperature fluctuations, determined by dedicated in situ 
measurements, are of the order of 10 _1 K/vTfz at I mHz, and dropping at higher 
frequencies within the MBW. The required stability conditions at the sensors, attached 
at the core's surface, thus need an attenuation factor of 1CP 5 , or better. Our 
analysis determines that a central aluminium core of 13 cm of radius, surrounded by 
a concentric layer of polyurethane 15-20 cm thick, comfortably provides the needed 
thermal screening which guarantees a meaningful test of the sensors' performance. 

The results of this paper are based on modelling. Because our aim is to produce 
a very stable thermal environment for the temperature sensors, we cannot check 
experimentally the correctness of our conclusions. We must instead rely on the validity 
of the hypotheses made — essentially that heat only flows by thermal conduction — and 
on the underlying physical laws which govern heat conduction. Even though there is 
good reason to believe that both are sufficiently accurate, unexpected behaviour e.g. 
at the interface between the metal core and the insulator, may partly distort the 
results. Direct measurements with very large temperature gradients applied across 
the insulating device are envisaged, and will be reported elsewhere as an auxiliary 
independent test of the model. 
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Appendix A. Thermal insulator frequency response functions 

Here we present some mathematical details of the solution to the Fourier problem, 
equations We first of all Fourier transform equations © and JJJJ: 

iuj pc p T(x, ui) — V • ^kVT(x,o;) (A.l) 

oo I 

T (9,<p;lj)=Y1 b lm (Lu)Y lm (0^) (A.2) 

i=0 m=-l 

Equation (|A.1(1 can be recast in the form 

(V 2 + 7l 2 ) f (x, u) = , < r < ai (A.3) 

(V 2 + 7 |) f (x, w) = , ai < r < a 2 (A.4) 
where r = |x|, and 

7 a = -«w i— , 7 2 = -iw ^— (A.5) 

Kl K 2 

To these, matching conditions at the interface" 1 " and boundary conditions must 
be added: 

f(r = a 1 -0,u))=f(r = ai+0,w) (A.6) 

df df 
Ki—(r = ai- 0, = k 2 — (r = ai + 0, u) (A.7) 
or or 

f(r = O2,u)=f (9,(p;u) (A.8) 

Equations (|A.3|I and (|A.4ll are of the Helmholtz kind. Their solutions are thus 
respectively given by 

{^2 Aim(w)3i('Yir)Y lm (8,(p) , < r < ai 

lm ( A -9) 
^ [Ci m (uj) jj(72r) + AmM ^(72^) ] *im(0, ¥>) , «i < i" < a 2 

where j; and yi are spherical Bessel functions 8 , 

, / 1 d V sin z , / 1 rf \ cos z ,. ... 

Ji(«) = 2 —j- — . yi{z) = -z l [ — - r ) (A.10) 

\ z dz ) z \ z dz ) z 

and the coefficients Ai m (u>), Ci m (u>) and Di m (u>) are to be determined by 
equations l|A.6l) ~ (|A.8() . These can be expanded as follows, respectively: 

J~] Ai m ji (7101) Y lm (6, <p) = 
= l C lm( w ) Ml2 a i) + Di m (u)yi(j 2 ai)] Y lrn {9, ip) (A. 11) 

lm 

Kiji ^ Ai m (u)fi(jiai)Yi m (0, ip) = 

lm 

= K272 l C i™(u) j'ti^ai) + Am(w)2/;'(72ai)] T Zm (<9,<^) (A.12) 
lm 

+ The temperature and the heat flux should be continuous across the interface. 



Theoretical foundations for. 



11 



^ [Cim(uj) .7/(7202) + Am (7202)] y/m(0,y) = 
£771 

= J] 6 im (^)y im (0,^) (A.13) 

Because of the completeness property of the spherical harmonics, the above 
equations completely determine the coefficients Ai m (uj), Ci m {u) and D/ m (w). The 
result is 

Ai m (uj) = £i(u>) h m (uj) , Q m (uj) = rn(u) 6/ TO (u>) , Am(w) = GM hm(w) (A. 14) 
with 

= a 7 \ [«272 j;(72ai)y/'(72ai) - k 2 72 i/'(72«i) 2//(720i)] (A.15) 

= a 7 \ [ K2 ^ 2 ii(7i«i)y/(720i) - K171 j/'(7iai) y/(72«i)] (A.16) 

= A 7 ^ I K i7i Ji(72ai) j;'(7iai) ~ K 272 j[ (7201) i/(7i«i)] ( A - 17 ) 
A/(w) 

and 

A/(w) = Ki7U'/(7i a i) [ii (7201)^/(7202) - i/ (72^2)2// (7201)] + 

+ «272i/(7i«i) [7/(7202) 2//(720i) - j/ (7201)^/(7202)] (A. 18) 

When the above results, equations l|A.15|) through (|A.18I) . are inserted back into 
equation (|A.9|) the result stated in equation 11U|) in the main text obtains, i.e., 

Hi m (x, u) bi m (uj) (A.19) 



where 



b(u)M'fir)Y lm (6,(p) , 0<r<ai 
H lm {x,co) = { (A.20) 
[^/(w)j/(72r) + 0(^)^(72^)] Y5 m (0,<p) , a 1 <r <a 2 

For monopole only boundary conditions, equation l|16|) . the transfer function is 

£oM jo(7i r ) , < r < 01 

H(r,w)=^ (A.21) 
i7oMio(72r) + Co M 2/0(72^) , 01 < r < a 2 
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